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Table of Scattering Function of Infrared Radiation 


for Water Clouds 


GUCHI YAMAMOTO, MASAYUKI TANAKA, AND SHOJI ASANO 


Geophysical Institute, Tohoku University, Sendai, Japan! 


ABSTRACT. Scattering of infrared radiation by water clouds has been calculated from 
the Mie theory. Tables are presented for: the phase function and the degree of polarization 
for each degree of scattering angle; the volume cross sections for extinction, scattering, 
and absorption; and the albedo for single scattering. Three models are used, representing 
altostratus, nimbostratus, and stratocumulus clouds. 


1. INTRODUCTION 


In the paper, “Radiative Transfer in Water Clouds 
in the Infrared Region” (Yamamoto, Tanaka, and 
Asano 1970), space did not allow the inclusion of the 
phase function for the clouds. Requests for these 
values have prompted us to prepare the tables 
published here. To satisfy general needs, the phase 
function is presented not only for the altostratus 
clouds discussed in our earlier paper, but also for 
nimbostratus and stratocumulus clouds. Volume 
cross sections for extinction, scattering, and 
absorption, and albedos for single scattering are also 
included. The theoretical basis for computation 
of the values is described below. 


2. SCATTERING THEORY OF LIGHT APPLIED 
TO CLOUDS 


2.1 Light Scattering by a Sphere 


Scattering of electromagnetic waves by a homo- 
geneous spherical particle is described by the Mie 
theory, and is discussed in detail by van de Hulst 
(1957) and by Kerker (1969). A summary of the 
results is shown below. 

The intensity functions i; and iz, which are the 
components perpendicular and _ parallel, respec- 
tively, to the reference plane (i.e., scattering plane), 
are given as follows: 


= 2n+] 
>» ere {an(x, m)77n(cos 0) 


n=1 


is(x,m,0) =| 


+ bu(x.m)ta(cos 0)) | (1) 


' Prepared under contract No. E-24-69(N). 


io(x,m,0) = S area {bn(x, m)Tn(cos 0) 


n=1 
+ an(x,m)Tn(cos ayy] (2) 


where an(x, m) and b,(x, m) are the Mie coefh- 
cients, and 77,(cos 8) and 7,(cos 8) are the angular 
coefhcients as defined by van de Hulst. The efh- 
ciency factors for extinction and scattering are 
given by 


Qext(x,m) = Si (2n+1)Re {an(x, m) 


+ br(x, m) } (3) 


Qsea(x, m) = 2 S (2n+1){ | an(x, m) i 


+ | ba(x, m) |?} (4) 


and the efficiency factor for absorption follows by 
subtraction, 


Oss = Oexi sz Osea : (5) 


Thus, the physical parameters involved are the 
scattering angle, 0, the complex refractive index of 
the sphere relative to the surrounding medium, m, 
and the size parameter 


x=— (6) 


where r is the radius of scattering particle, and 
d is the wavelength of light in the surrounding 
medium. 


2.2 Light Scattering by Clouds 


The total scattering and absorption cross sections 
per unit volume of the clouds are obtained by inte- 
grating the Mie single-particle cross section with 
respect to the appropriate size distribution. The 
volume cross section, o@, may be given by the 
integral 


o(m,A) = ie ar7Q Ge m) n(r)dr (7) 


where n(r)dr is the number of particles per unit 
volume with radii between r and r+ dr, Q is Qezt, 
Qsca, OF Qavs, and r; and rz are the lower and upper 


limits, respectively, to the radii of the particles. 


The albedo for single scattering, w, is defined by 


Osca 
(Osca + ads) 


__ Osca 


@ = 


(8) 


Oext 


where Gsca, Gads, and Gert are the volume cross 


sections, respectively, for scattering, absorption, 
and extinction obtained from (7). 


The elements of the phase function normalized 
to unity are given by the integrals 


P;(m, 0, d) = ( A i Ji im n(r)ij(r, m, 0)dr 


] 
4ar 27 Osca Jry 


G=1,,2)~ (9) 


where j= 1 and ;=2 correspond to the perpen- 
dicular and_ parallel components, respectively. 
The phase functions are normalized to satisfy the 
condition: 


| P(cos 0) oe ] (10) 


where P(cos @) is the averaged phase function given 
by 
P(cos 6) = 2[P,(0) + P2(6)]. (11) 


For unpolarized incident radiation, the scattered 
radiation shows partial linear polarization, with the 
maximum electric field oscillation being either 
perpendicular or parallel to the scattering plane. 
The degree of polarization is given by the ratio 


P,(0) race P.(0) 


P,(6) + P2() 2 


with positive and negative values corresponding, 
respectively, to the two cases mentioned above. 


3. TABULATED FUNCTIONS 


The optical properties of water clouds are shown 
in tabular and graphical forms in the wavelength 
region from 5 to 40 wm. In the calculation, we used 
the value of complex refractive index of liquid water 
evaluated by Pontier and Dechambenoy (1966). The 
drop size distributions of our cloud models (fig. 1) 
are those of altostratus (As), nimbostratus (Ns), and 
stratocumulus (Sc) measured by Diem (1948). 

The intensity functions (1) and (2) and the efh- 
ciency factors for extinction (3) and scattering (4) 
were computed using the computational scheme and 
recursion criterion of Deirmendjian, et al. (1961). 

To evaluate the integrals (7) and (9), we replaced 
each size distribution curve shown in figure 1 by 
histograms, using radius increments, 6r, of 0.25 um 
(Carrier et al. 1967); particle concentrations for 
each radius, rx, are given by n(rx) drx. The validity 
of our calculation was checked by comparisons 
with the published results of Deirmendjian (1964), 
Bauer (1964), and Carrier et al. (1967). 

The symbols used in the tables are as follows: 


EXT Volume extinction cross section, 
Oext.incm!. 

SCA Volume scattering cross section, 
Osca,incm=}, 

ABS Volume absorption cross section, 
Gaps, i CMa", 

ALBEDO Albedo for single scattering, w, per 
unit volume. 

ANGLE © Scattering angle, 0, in steradians 
and in degrees measured from 
the direction of propagation of 
incident light. 

PHASE Normalized phase function, P(cos 8), 
per unit volume, defined by (11). 

POL Degree of polarization, defined by 


(2). 


4. SOME OPTICAL PROPERTIES OF WATER 
CLOUDS 


In the infrared region from 5 to 40um, liquid 
water has absorption bands (Yamamoto et al. 1970) 
near 64m (the v2 band) and 15um; the latter is a 
broad, strong band caused by intermolecular 
vibration, which does not exist in the vapor phase. 

Figure 2 shows the spectral distribution of volume 
extinction, scattering, and 


cross. sections for 
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FicurE 1.—Droplet size distribution of various cloud models. As, Ns, and Sc denote altostratus, nimbostratus, and stratocumulus, 
respectively. The number density and liquid water content are also shown in the figure. 


absorption of the altostratus model. The absorption 
cross section is large in the band regions, whereas 
the scattering cross section depends mainly on the 
real part of the complex refractive index. At wave- 
lengths shorter than 10.5um, the scattering cross 
section is larger than the absorption cross section; 
in the wavelength region beyond 10.5um, the re- 
lationship is reversed. The total extinction cross 
section is the sum of the scattering and absorption 
cross sections. 

Figure 3 shows the spectral distribution of volume 
extinction cross sections (solid lines) and of albedos 
for single scattering (broken lines) for the three 
nimbostratus, and 


cloud models  (altostratus, 


stratocumulus). The extinction cross section is 
largest in nimbostratus, which, as shown in figure 1, 
contains larger cloud droplets than the other clouds 
and is optically the most dense (highest liquid water 
content); it is smallest in the optically thin (lowest 
liquid water content) stratocumulus. In the far 


infrared region beyond about 10.5um, albedos of 


the three cloud types are less than 0.5; that is, the 
contribution of absorption to the total extinction 
is larger than that of scattering. 

Figure 4 shows the values of the phase function, 
Pi(cos 0), for the altostratus model at wavelengths 
of 5, 10, 20, and 40um as a function of scattering 
angle, 0. The predominance of the forward scatter- 
ing, which is a characteristic property of Mie 
scattering, is clearly exhibited at all wavelengths. 
However, it is more marked at shorter wavelengths, 
for which the size parameter x (i.e., particle size 
relative to wavelength) becomes larger for a given 
size distribution of cloud droplets. The pattern of 
the phase function would become flatter and 
smoother in a cloud composed of droplets with 
smaller size parameters. 

These features also are evident in the comparative 
phase functions of altostratus, nimbostratus, and 
stratocumulus at a wavelength of 10mm (fig. 5). The 
forward scattering is most marked in nimbostratus, 
which is composed of droplets with the largest 
size parameters among the three model clouds. 
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FIGURE 2.— Volume extinction (o,,.1), scattering (Gyeq), and absorption (y»s) cross sections of altostratus as a function of wavelength. 
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FIGURE 3.—Volume extinction cross section (solid lines) and the albedo for single scattering (broken lines) for nimbostratus (Ns), alto- 
stratus (As), and stratocumulus (Sc) as a function of wavelength. 
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FIGURE 4.—Phase function P(cos 6) for altostratus at wavelengths of 5, 10, 20, and 40mm. 
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FIGURE 5.— Phase function P(cos @) for stratocumulus (Sc), altostratus (As), and nimbostratus (Ns) at wavelength of 10 wm. 
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